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We study the zero-temperature (T = 0) quantum rotor model with on-site disorder in the charging 
energy. Such a model may serve as an idealized Hamiltonian for an array of Josephson-coupled 
small superconducting grains, or superfluid 4 He in a disordered environment. In the approximation 
of small-amplitude phase fluctuations, the Hamiltonian maps onto a system of coupled harmonic 
oscillators with on-site disorder. We study the effects of disorder in this harmonic regime, using 
the coherent potential approximation (CPA), obtaining the density of states and the lifetimes of 
the spin- wave-like excitations for several choices of the parameters which characterize the disorder. 
Finally, we estimate the parameters characterizing the T = quantum melting of the phase order, 
using a suitable Lindemann criterion. 
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I. INTRODUCTION 

The quantum rotor model is a widely studied Hamil- 
tonian which may serve as an idealized representation 
of many physical systems. It consists of two parts: a 
potential energy which represents the coupling between 
two-component spin variables of fixed length (an "XY" 
Hamiltonian) and the kinetic energy of those spins. The 
model is quantum-mechanical because the potential and 
kinetic energies do not commute. 

One possible realization of this model is a Josephson 
junction array, i.e., a collection of superconducting is- 
lands with Josephson or proximity effect tunneling be- 
tween the islands. Such arrays can be produced experi- 
mentally in a wide range of geometries, and with a great 
variety of individual junction parameters. They may be 
described by the quantum rotor model in the limit of 
small grains. In this case the potential energy corre- 
sponds to the Josephson coupling, and the kinetic en- 
ergy corresponds to the charging energy. Such arrays 
can serve as model systems for studying quantum phase 
transitions^ under conditions such-,that the experimental 
parameters can be readily tuned.cl Recently, it has also 
been proposed that Josephson junctions may serve as a 
quantum bit (qubit) in quantum information technology, 
where the quantum logic operations would be performed 
by controlling gate voltages or magnetic nelds.0 

The quantum rotor model may also be realized in other 
physical systems. For example, in liquid 4 He, the poten- 
tial energy may represent a discretized form of the gra- 
dient energy term in the Ginzburg-Landau free energy of 
superfluid 4 Hc, and favors a state in which all parts of 
the superfluid have the same phase, while the kinetic en- 
ergy term is the energy associated with the time- variation 
of the order parameter phase. The model may be par- 
ticularly useful for helium films and helium in confined 
geometries such as porous glasses. The spectrum of ele- 
mentary excitations in such confined geometries has re- 
cently been studied extensively.! 4 ! This spectrum may be 
treatable by the quantum rotor model. 

In this paper, we study a special kind of quantum ro- 



tor model, namely, the two-dimensional model with di- 
agonal disorder in the kinetic energy. Our emphasis in 
this study is to understand what happens to the spec- 
trum of elementary excitations in the presence of this 
disorder. We also investigate the conditions for the de- 
struction of phase order in this model arising from such 
quantum fluctuations in the phase variables. This model 
includes a particularly simple form of disorder. In the 
case of a Josephson junction array (JJA), it would corre- 
spond to fluctuations in the self-capacitance of individ- 
ual grains, but no significant randomness in the Joseph- 
son coupling between the superconducting grains. In the 
case of superfluid helium, the diagonal disorder might 
correspond to randomness in the pore geometry (for a 
three-dimensional system) or arising from surface rough- 
ness (for a thin film). 

In the present paper, we are primarily interested in 
the spectrum of elementary excitations in this model at 
low temperatures. We therefore adopt a simple mean- 
field approximation which has been successfully used to 
treat elementary excitations in other systems. This is the 
coherent potential approximation|-.£CPA), as described 
originally by Soven,El Taylor et alfm and Velicky et al.tt 
For the present problem, the CPA may be adequate at 
very low temperatures, under conditions such that the 
original Hamiltonian can be approximated as harmonic. 
At higher temperatures, or under conditions such that 
the harmonic approximation is inadequate, it will be 
necessary to consider a more sophisticated approxima- 



tion, such. 
(SCHA),!' 1 



Jie self-consistent harmonic approximation 
or a field-theoretic treatment in which the 
charging energy is mapped onto a coupling in the imagi- 
nary time direction.E£l 

The remainder of this paper is organized as follows. In 
the next section, we present the model Hamiltonian, and 
its relation to a set of coupled harmonic oscillators. In 
Sec. Ill, we describe how the CPA may be used to study 
the disorder-averaged properties of the model Hamilto- 
nian in the harmonic approximation. Our numerical re- 
sults based on the CPA are presented in Sec. IV, followed 
by a concluding discussion in Sec. V. 
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II. MODEL HAMILTONIAN 

We will consider the following model Hamiltonian: 

H = \ E u i < + E - cos ^ - w 

i jk 

Here 9j is the phase of the quantum rotor, and rij is a 
variable which is quantum-mechanically conjugate to the 
angle 8j. rij and 9k are assumed to have commutation 
relations 



-id 



jki 



(2) 



where Sjk is a Kronecker delta function and [. . .] repre- 
sents a commutator. In Eq. ([!]), the first sum runs over 
sites j, and the second sum runs over all distinct pairs 
of sites j and k. If this model is assumed to describe a 
Josephson junction array, then the first term represents 
the capacitive energy of the array in the diagonal ap- 
proximation. The most general form of this capacitive 
energy would be | J2jk ^jk n j n k't m Eq- ([j]), we make the 
approximation Ujk = UjSjk, where Uj is related to the 
diagonal capacitance Cj, or alternatively to an effective 



"mass" M. 
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The second term represents the sum of the Josephson 
coupling energies between grains j and k. In Eq. (^), 
it is implicitly assumed that there is no contribution to 
H arising from dissipation. In this paper, we shall as- 
sume that Jij ^ only if i and j are nearest-neighbor 
grains, and that, furthermore, there is no disorder in this 
nearest- neighbor coupling, i.e., that Jij — J for all near- 
est neighbor pairs i and j. Hereafter, we shall refer to 
this model Hamiltonian as describing a Josephson junc- 
tion array, although, as noted above, the model is also 
applicable, in principle, to other physical systems. 

For sufficiently low temperatures, we assume that the 
above Hamiltonian can be approximated by its harmonic 
form: 



larm — ^ 
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9 k ) 2 =K + V, (4) 



where the second sum runs over distinct nearest neighbor 
pairs. While one can go beyond the harmonic approxima=j 
tion by using a self-consistent phase phonon approach,Bt£l 
we shall here consider only the harmonic approximation, 
which is more easily combined with a treatment of dis- 
order. The Heisenberg equation of motion of the phase 
variable is iTiOj = [6j,H] = iUjUj, where the last equal- 
ity follows from the commutation relation (^). This equa- 
tion makes clear that -f/harm is formally analogous to an 
array of "masses" Mj harmonically coupled by springs of 
"spring constant" J. However, in contrast to a real lat- 
tice of masses and springs, there is no polarization degree 



of freedom. K and V represent the kinetic and potential 
energy terms in this Hamiltonian. A disordered array of 
such oscillators can thus be treated by the same methods 
used to treat disordered arrays of coupled oscillators. 

In this paper, we shall consider a bimodal distribution 
of the Uj 's. Thus we take Uj — U\ or U2 with probabil- 
ities p and 1 — p respectively. We also assume that the 
type 1 and type 2 grains are randomly distributed on an 
ideal lattice, with no correlation between the occupan- 
cies of neighboring sites. With this choice, the Hamilto- 
nian (Q) is analogous to the classic problem of substitu- 
tional mass disorder in phonon spectra, except that the 
"phase phonon" excitations have only one possible po- 
larization. This problem has previously-.been studied by 
Sovenjj Taylom, and Elliott andflay lor ,Q using the coher- 
ent potential approximation.BEtLlu In the present paper, 
we shall use the same approach to study the effects of 
disorder in the "phase phonon" spectrum. Note that be- 
cause the entire Hamiltonian is invariant under a global 
rotation of phases, the long-wavelength phase phonons 
are Goldstone bosons with frequencies that go to zero 
linearly in the wave vector k, even in the presence of 
disorder, and even in two dimensions. 



III. COHERENT POTENTIAL 
APPROXIMATION 

A. General equations 

The CPA provides an approximate way of calculating 
the disorder-averaged Green's function 



G(z) 



1 



harm / dis 



(5) 



where z is a complex variable, -f/harm is the Hamiltonian 
operator, and (...}dis represent an average over disorder 
realizations. Of particular interest is the diagonal matrix 
element of this operator, which we write 



F{z) = {3\G{z)\j) = G{j,r,z). 



(6) 



Because the operator G is disorder-averaged, F(z) is in- 
dependent of j. 

Also, since G(z) is disorder-averaged, it can be ex- 
pressed in terms of a periodic operator H c r by the re- 
lation 



G(z) 



1 



z-H, 



cir 



(7) 



For the present case of "site-diagonal" disorder, H e s can 
be written 



where K e g can be expressed as 



2M eG ^ J 



(8) 



(9) 
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Here M c g(z) is a complex, frequency-dependent "effec- 
tive mass" which is, however, the same for each site j. 
It is convenient to express M e g in terms of two other 
quantities - the "virtual-crystal" mass 



M 



VGA 



P M 1 + (l-p)M 2 , 



(10) 



and a dimensionless, complex self-energy function E(z) 
defined by 



M off (z) = Mvca [1 - E(z)] . 



(11) 



In the CPA, it is found that F(z 2 ) pad E(z) are con- 
nected by the complex scalar equatior£3 



MvcaS(z) 



p(l -p)(SM) 2 z 2 F(z 2 ) 



1 + z 2 [(1 - 2p)SM + MvcaS(z)] F(z 2 ) ' 

(12) 

Here SM = M 1 - M 2 = -e M 2 , where e = 1 - M 1 /M 2 = 
1 — C1/C2, and we have adopted the convention that 
Ci < G 2 . 

In order to complete the approximation, we need to 
express F(z 2 ) in terms of E(z). This is readily ac- 
complished as follows. First, we consider the analo- 
gous virtual-crystal Green's function Gvca(z) = l/(z — 
■frVcA - V), where K VC A = J2 3 - n 2 /[2M VC a]- Since 
KycA + V is periodic, Gvca(z 2 ) is diagonal in momen- 
tum space. In particular, the diagonal matrix element is 
given by 



(q|G VC A(^)|q) 



1 



= G V CA(q;2 2 ), (13) 



where cj q is the frequency of a "phase phonon" of wave 
vector q in the virtual crystal, and |q) and \i) are related 
by |q) = N~ 1 / 2 X)j e<q ' Ri K)> being the position of 
the ith grain and N being the number of grains in the 
lattice (which may be assumed to have periodic boundary 
conditions) . The corresponding diagonal matrix element 
Fvca(z 2 ) is given by 



F VC a(z 2 ) - 



1 



2NM VCA 



^G VCA (q,z 2 ), (14) 



where the sum runs over q's in the first Brillouin zone of 
the grain lattice. Next, we introduce the virtual-crystal 
density of states <?vca(?7) by 



9vca{v) = S (V ~ 



in terms of which Fvca(z 2 ) takes the form 



(15) 



■-VCA 



1 



Mvca 



9vca(v), 

— 2 dr >- 



Finally, we can obtain F(z 2 ) for the actual disordered 
crystal in terms of Fvqa{z 2 ) simply by making the re- 
placement z 2 — > z 2 [1 — S(z)], so that 

F(z 2 ) = F VCA (z 2 [1 - E(*)]) 

1 f°° 9vca(v) 



Mvca J z 2 [1 - E(*)] - V 



■drj. (17) 



Equations. ( pL 2| ) and (17) represent two coupled equations 
for F(z 2 ) and E(z) which can be solved simultaneously, 
given a model for svca(^). 

In the present paper, we will consider a two- 
dimensional (2D) lattice with a simple Debye-like density 
of states defined by 



9VCa(v) 



'/ 



2-kv 2 



@(u>D - »/)• 



(18) 



Here v = \J J/A/vca is the speed of the "phase phonon" 
excitations in the virtual crystal, and Q(ll>d — ?y) is a unit 
step function. The Debye cutoff frequency 



UJ D 



(47TW 2 ) 1 / 2 = (47tJ/A/ V ca) 1/2 ; 



(19) 



follows from requiring that the number of phase phonon 
states be normalized to one per grai n. The quantity 
Fvca(z 2 ) follows readily from Eq. (jig): 



Fy CA (z 2 



1 



In 



■ UJ, 



(20) 



MvCAwfj 

Equations (|l2|), (Jl7|) , and ( p0| ) constitute a self-consistent 
set of equations which can now be explicitly solved to give 
the properties of the phase excitations in the disordered 
system, within the CPA. 



B. Lattice properties 

1. Density of states 

The density oL phase phonon excitations within the 
CPA is given bytL^I 



ScpaM 



■Im 



V 2 9vca(v) 
; 2 [1 -E(w)] -v 2 



dr]. (21) 



Using Eq. (|1 
simplified to 



for gvCA(f?), this expression can be further 



ScpaM = lm[h(a 



TTUl 



(22) 



where x 2 
x 2 In [x 2 /(x 2 



= (^M) [1 
- 1)] - 1. Thus, 



E(w)] and h(x 2 ) = 
once we have obtained 
the self-energy from Eq. (|12|), we can find the density 
of vibrational states, <?cpa(w) from Eq. (p2|).E-3 

The CPA can also be used to give the component den- 
sities of states g^?p A (uj) (a = 1, 2), defined to be the aver- 
age density of phase phonon density of states on grains of 
type a. This is obtained in terms of the average diagonal 
(16) matrix element (j, j; z 2 ) for a grain oLtype a. In the 



CPA, this matrix element takes the formfl 



F(z 2 ) 



1 - {Ah 



VCA 



M a -M VC A^)z 2 F(z 2 Y 
(23) 

where F(z 2 ) is defined in Eq. (§). The average density of 
states at a grain of type a (normalized so as to integrate 
to unity over positive frequencies) is given by 
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( Q ) i \ 



-Im 



h^ 2 ) 

1 - (Mvca -M a - M VCA E)cu 2 F(uj 2 ) 



(24) 



where cc 2 and ft(ic 2 ) were defined below Eq. (j22j). 

Within the CPA, the component densities of states have the pleasing property 

PScpaM + ( 1_ TOa( w ) = 5cpa(w). 

Thus, the total phase phonon density of states equals the sum of the properly weighted partial densities of states 
on the two species of grains. In our calculations, we have confirmed that this property is, indeed, satisfied for our 
particular choice of the Debye model. 



(25) 



2. Spectral Function 

Our self-consistent set of equations permits calculation of another useful quantity: the spectral function, defined by 



ImG(q,0 2 ) = 



'■Im [S(w 2 )] 



{u 2 {\ - Rc [£(^ 2 )]) - V} 2 + {w 2 Im [E(w 2 )]} 2 ' 



(26) 



The spectral function gives the frequency distribution of 
excitations of wave number q in the array, and, for our 
Debye model, depends only on the magnitude q. The full 
width at half maximum (FWHM) of the spectral function 
is inversely proportional to the decay time of the spin- 
wave-like excitations at wave vector q. 



3. Mean-Square Phase Fluctuations 

Finally, we can use our CPA calculations to infer 
the mean-square phase fluctuations at the jth grain. 
We denote this quantity {\9j(t)\ 2 ), where (...) denotes 
a quantum-mechanical average, and 9j(t) denotes the 
phase of the jth grain at time t. A quantum-mechanical 
average is appropriate at temperature T = 0. At finite 
temperatures, one should carry out both a quantum and 
a thermal average. However, such an average diverges at 
T ^ in two dimensions (d = 2), because of contribu- 
tions from the long- wavelength phase phonons. In the 
present discussion, therefore, we discuss only the T = 
limit. 

We will consider the quantity (\9j(t)\ 2 ) for a lat- 
tice with Josephson coupling J and self-capacitance Cj, 
within the harmonic approximation. In this case, one 
can simply adapt the discussion of Ref. |l5| for lattice vi- 
brations in mass-disordered systems. In fact, the present 
problem is slightly easier since there are no polarization 
degrees of freedom. To obtain this quantity, we first use 



a result from Ref. h 5 that 



;(t)0 fe (O))=--Im / 

71" J-, 



'GO',*^ 2 )*, (27) 



where, at T = 0, G(j,fc;0 2 ) = The desired 

quantity is actually \\m t ->o(8j(t)9j(0)) averaged over dis- 
order realizations, which we write 



9? I )dis 



— Im / F{0 2 )d0, 
t Jo 



(28) 



where we have used the definition of F{z 2 ) in Eq. (||). 
The operators 0j are computed at times t = (though 
this average is time- independent). 

We first consider the virtual crystal lattice, in which 
Cj = Gvca for all j . In this case, (|%| 2 ) is readily calcu- 
lated from Eq. (j28|), using Fvca(z) from Eq. (|2(]|), with 
the result 



h 



/VCA 



(29) 



yJl6we 2 J/{h 2 CvcA). 

we can calculate 



where 0d = y^rJ/MvCA 

For the actual disordered lattice, 
within the CPA not only the full disorder average of 
(|#j| 2 )dis, but also the disorder average of (|#j| 2 ) over 
sites of type a, which we denote {\9j\ 2 )^ a K The result is 



9 |2 } («) = _l lm j G<*(j,r,0 2 )d0 C'»0! 

7T 



-h f°° F(x 2 ) 

7r Im J 1 - (Mvca - Af Q - M VCA S)0 2 ^ 2 ) ' ( } 
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FIG. 1: (a) The real (dotted line) and imaginary (dashed 
line) parts of the self-energy, £(0), plotted as a function of the 
scaled frequency 0/0d for three different cases: e = 0.1,0.5, 
and 0.9, at a concentration p = 0.25 of light defects. Note 
that ujd is a function of p and e [cf. Eq. (J19|) ] - (b) Same as 
(a) except that p = 0.75. Note the different vertical scales on 
different portions of the Figure. 



where x 2 = [uJ 2 juJ 2 ^) [1 - 
the CPA satisfies 

(l^| 2 )d is (p) 



E(w)]. The full (|^| 2 )dis(p) in 



(32) 



where the sum runs over the two species a, and in our 
notation p\ = p, and pi = 1 — p. 

We can use these equations to compute various ratios, 
such as (|#j| 2 )dis(p)/(|#j| 2 )vCA, and hence to see how the 
mean-square phase fluctuations increase with increasing 
disorder. 



IV. NUMERICAL RESULTS 

We have solved the model described in the pre- 
ceding section to obtain the self-energy, density of 
states, spectral function, and mean-square phase fluctu- 
ations of a two-dimensional lattice, using a Debye phase 
phonon spectrum. The numerical solution are obtained 
from straightforward iterative solution of the three self- 
consistent equations (fl2]), (p7|), and (|20|). It is important 
to take care that the resulting solutions are the proper 
physical ones, with ImS(w) < and with both F(z 2 ) 
and E(z) varying smoothly with z. 
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FIG. 2: (a) Normalized density of states of the phase phonon 
excitations, plotted as a function of the scaled frequency 0/0d 
calculated for a concentration p = 0.25 of low-capacitance 
defects and three different values of the disorder parameter 
e = 1 — C1/C2, as indicated in the Figure. The total density of 
states for phase phonon excitations, (;cpa(0), is the solid line, 

while the component density of states <7cpa(0) ano - ScpaM 
associated with low-capacitance and high-capacitance grains 
are the dashed and dotted lines, respectively. The concen- 
tration of light defects (type 1) is p = 0.25. (b) same as (a) 
except that p — 0.75. As in Fig.l, u)d is a function of both p 
and e [Eq. (Ra)]. 



In Fig. 1 we show the real and imaginary parts of 
the self-energy, £(0) (dotted and dashed curves respec- 
tively), plotted against the scaled frequency 0/0_d, where 
0d = (47tJ/Mvca) 1 ^ 2 is the virtual crystal Debye fre- 
quency. We have carried out our calculations for two 
concentrations of the low-mass defects (type 1), p = 0.25 
and p = 0.75; the results are shown in Figs. 1(a) and 
1(b). The three different panels correspond to different 
values of the ratio e = 1 — Ci/C 2 = 0.1,0.5, and 0.9. 
These ratios range from small deviations in capacitance 
between type 1 and type 2 (e = 0.1) to large deviations 
(e = 0.9). Note the differences in the scales of the plots. 

Given the self-energy, several other physical quantities 
can be calculated, as explained in Sec. III. For example, 
the averaged density of states of phase-phonon excita- 
tions can be obtained within the CPA from Eq. (|2^), 
while the averaged component densities of states follow 
from Eq. (p4|). In Fig. 2 we show the average density of 
states gcPAX 10 ) obtained within the CPA (solid line), as 
well as the component densities of states 5cpa(0) (dashed 
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FIG. 3: (a) Spectral function ImG(q, 2 ), plotted against the 
scaled frequency 0/0d for three different mode frequencies 
uiq as indicated in each panel. The concentration of light 
defects is p = 0.25. For each tj q we have considered e = 
1-Ci/Ca =0.1 (solid line), 0.5 (dotted line), and 0.9 (dashed 
line), corresponding to increasing amounts of disorder. For 
0q/0D < 1, the spectral function always exhibits a single 
sharp peaked centered at = q . (b) same as (a) except that 
p = 0.75. 



( 2) 

line) and <?cpa(0) (dotted line). The sum rule given in 
Eq. (|25|) is evidently satisfied by these three densities of 
states. 

These results show characteristic features expected 
from CPA calculations. For the case where the capac- 
itance of type 2 is close to that of type 1 (e = 0.1), the 
full and the two component densities of states are all very 
similar, and all resemble the virtual-crystal result. As e 
increases, there start to be more conspicuous differences 
between the partial densities of states on the light-mass 
(type 1) and heavy-mass (type 2) grains. In particular, 
there is a clearly-developed band gap between the two 
classes of states for e = 0.9. This behavior is well known 
as the "split-band" regime in the phonon problem. The 
light-mass phase phonons presumably correspond to lo- 
calized modes, though this localization is not probed in 
the CPA. 

In Fig. 3, we plot the spectral function defined in Eq. 
( p6| ) as a function of the scaled frequency 0/0d for several 
choices of parameters. Our calculations are carried out 
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FIG. 4: The full width at half maximum (FWHM) of the 
spectral function, plotted versus scaled frequency 0/0d for e = 
0.2, and for several concentrations p of the small-capacitance 
grains, as indicated in the legend. 



at three different mode frequencies, q /0D = 0.2, 0.8, 
and 1.2, which are shown in the three different panels, 
and for two different defect concentrations, p — 0.25 and 
p = 0.75, which are shown in Figs. 3(a) and 3(b). In 
each of these cases, we have also considered three val- 
ues of the mass ratio parameter: e = 0.1,0.5, and 0.9. 
(Note in particular the vertical scale changes in the Fig- 
ures.) For q /0D < 1, the spectral functions in all cases 
exhibit single sharp delta-function-like peaks centered at 
= q . However, for q /0_o > 1, there is some spectral 
strength which shows up as weak peaks at higher frequen- 
cies; these are presumably due to the localized modes in 
the light-mass bands. Note also that for e = 0.1 there is 
almost no spectral weight for q /0n > 1. 

In each of the cases of Fig. 3, we have also calcu- 
lated the full width of the spectral function peak at half- 
maximum (FWHM) . This width is inversely proportional 
to the lifetime of the phase phonon like excitations. In 
Fig. 4, we plot this FWHM as a function of q /0_o for 
e = 0.2, and for several concentrations of the light defects 
as indicated in the legend. Evidently, for a given w q , the 
excitation lifetimes decrease, as expected, as p increases 
from to 0.5, then increases again as p increases from 
0.5 to 1. This behavior is consistent with the well-known 
result of second-order perturbation theory that this life- 
time should vary approximately as l/[p(l — p)] in the 
weak-disorder rcgime.u 

Finally, we turn to the transition from phase coherence 
to phase incoherence in this system. If the model Hamil- 
tonian is assumed to represent an array of small Joseph- 
son junctions, this corresponds to a superconductor- 
insulator (S/I) transition. It is believed that supercon- 
ductivity in Josephson junction arrays is destroyed when 
the ratio of the charging energy to the Josephson cou- 
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FIG. 5: (a). Mean-square phase fluctuation (\9j | 2 )dis/(|#j | 2 )vca, as evaluated in the CPA, plotted as a function of e = 1 — C1/C2 
for p = 0.1. The solid line shows the CPA average value, the dashed line shows the mean-square phase fluctuations on a type 
1 grain, and the dotted line on a type 2 grain, (b). same as (a) except that p — 0.3. (c). same as (a) except that p = 0.5. 



pling energy exceeds some characteristic limiting value.BI 
In fact, the S/I transition is a standard example of a 
T = quantum phase transition.!!! In two-dimensional 
square arrays with nearest-neighbor coupling, it has been 
found experimentally that the S /I transition occurs near 
U/J = 1.7 = (U/J) CI E 

In our discussion, we will assume a simple Lindemann 
melting criterion, namely, that superconductivity is de- 
stroyed when the mean fluctuations in the phase, given 
in Eq. exceeds some limit to be determined below. 

(In d = 2, such a criterion can be plausible only at T = 0, 
since the mean-square phase fluctuation diverges at any 
finite T '.) If all the grains have the same charging energy, 
then the mean-square fluctuation is given by Eq. (29), 
which implies a critical value of 



% = (U/J)H 2 « 1.3. 



(33) 



Here j can be any lattice site, and we have used 
(U/J) CI = 1.7. 

How can this Lindemann criterion be extended to a 
lattice with diagonal disorder in the charging energy? 
We speculate that a modified Lindemann criterion may 
still be usable, in the following way. Let us consider the 
model studied in this paper, in which there are two types 
of grains: "light" and "heavy" (or small and large capac- 
itance). We suggest that a possible Lindemann criterion 
for destruction of superconductivity in this case is that 
the mean-square phase fluctuations, averaged over the 
two types of grains, should exceed the same critical value 
as in the ordered case. 

In Fig. 5 we plot the calculated ratio 
{\0i\ 2 )^ /{\9i\ 2 )vcA, which gives the mean-square 
phase fluctuations on grains of type a — 1,2 with 
respect to those of the virtual crystal (dashed and 
dotted lines). We also show (|#i| 2 )dis/(|#i| 2 )vcA which 
represents the average for the entire lattice (solid line). 
All quantities are calculated in the CPA as a function 
of e = 1 — C1/C2 for several different values of p, as 
indicated in the Figure. As expected, and as is clear 



from the Figure, for all values of p, (|#i| 2 ) (1) > {\6i\ 2 )^ 
because component 1 represents the light defects. 

In principle, one might make an estimate of the criti- 
cal value of {U/J) 1 / 2 at which melting occurs, as a func- 
tion of p and e, by assuming that melting occurs when 
(\9i\ 2 ) > 1-3, as in the ordered lattice. Thus, given the 
properties of the lattice, one can easily extract an esti- 
mate of the melting parameters from this Figure. 



V. SUMMARY 

In this paper we have considered the behavior of a 2D 
quantum rotor model with diagonal disorder. In the case 
of Josephson junction arrays, this model corresponds to 
an array of small Josephson junctions with diagonal ca- 
pacitive disorder. We have used the CPA to estimate the 
effects of disorder on the phase phonon density of states, 
the spectral functions, and the self-energy, all within the 
harmonic approximation to the original Hamiltonian. Fi- 
nally, we obtained a crude estimate of the parameters 
governing the transition from coherence to incoherence 
in the disordered system, using a simple Lindemann cri- 
terion which may be applicable at T = 0. 

Viewed as a representation of a Josephson junction ar- 
ray (or a thin superconducting film), our model, and our 
approximations used to treat it, are both very simpli- 
fied representations of a real superconducting array. For 
example, a bimodal distribution of charging energies is 
oversimplified, as is our assumption of diagonal charg- 
ing energies and diagonal charging disorder. In addition, 
we have not considered dissipation, which is known to 
have important effects on S/I transitions in Josephson 
arrays.Ell Nevertheless, our results may provide some use- 
ful insights in understanding the S/I transition in such 
arrays, as well as the spectrum of excitations to be ex- 
pected in disordered arrays. 

Would it be possible to measure experimentally the 
phase phonon spectrum we have calculated in this pa- 
per? It is not clear what experiment would be directly 
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sensitive to this spectrum in a superconducting array or 
film. A more promising direction might be the excita- 
tion spectrum of a Bose superfluid (e.g. 4 He) in a porous 
medium. Such spectra have been extensively studied ex- 
perimentally r -jp.rim ari ]y using inelastic neutron scattering 
techniques .□llj'EZrHEj While the behavior of 4 He in such 
porous media is certainly more complex than the rela- 
tively simple model discussed here, some aspects of the 
observed behavior (e. g., the persistence of rather sharp 
excitation peaks even in highly disordered systems) seems 
to be mirrored in our calculated spectral functions. Ob- 
viously, a more refined model, aimed specifically at the 
geometries of these porous glasses, is needed before any 
comparison to experiment can be contemplated. 

Finally, we briefly comment on the quality of the ap- 
proximation itself. While the CPA is an excellent mean- 



field approach, it is still based on a harmonic approxi- 
mation to the underlying Hamiltonian. Ideally, it would 
be preferable to obtain the desired spectral functions us- 
ing a more accurate approach, such as a quantum Monte 
Carlo technique (see, e.g., Ref. [2(]), niL.a_version of the 
self-consistent phonon approximationjijiid suitably gen- 
eralized to treat a disordered system. We hope to return 
to such approaches in a future publication. 
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